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Abstract 

The recent paper by W.-X. Wang, Y.-C. Lai, J. Ren, B. Li & C. Grebogi |arXiv:1107.2T77v l] Q pro- 
posed a method for the control of complex networks with nonlinear dynamics based on linearizing 
the system around a finite number of local desired states. The authors purport that any bidirec- 
tional network with one-dimensional intrinsic node dynamics can be fully controlled by a single 
driver node, which can be any node in the network, regardless of the network topology. According 
to this result, a network with an arbitrarily large number of nodes (a million, a billion, or even 
more) could be controlled by a single node. Here we show, however, that this result is specious 
and that it does not hold true even for networks with only two nodes. We demonstrate that the 
erroneous results are a consequence of a fundamental flaw in the proposed method, namely, that 
reaching a local desired state using the proposed linearization procedure generally requires a non- 
local trajectory, which grossly violates the linear approximation. Their further conclusion that 
network systems with nonlinear dynamics are more controllable than those with linear dynamics 
is a known fact presented as novel based on a flawed argument. A central problem underlying 
the authors' argument is that their formulation is entirely based on the number of driver nodes 
required to reach a desired state — even when this number is one, keeping the system at that de- 
sired state generally requires directly controlling all nodes in the network. When this is taken into 
account, the conclusion that nonlinear dynamics can facilitate network control was already antici- 
pated in S.P. Cornelius, W.L. Kath & A.E. Motter |arXiv:1105.3726V l] 0, which is not referenced 
in their paper. If one insists on using the authors' formulation, then, in contradiction to their 
claim, nonlinear systems would never be more controllable than linear ones simply because the 
controllability of linear systems requires a single driver node, as shown by N.J. Cowan, E.J. Chas- 
tain, DA. Vilhena & C.T. Bergstrom |arXiv:1106.2573V l] 0, which is also not referenced in their 
paper. Taken together, all results in Ref. [l| are unsound, including the authors' generalizations to 
high-dimensional node dynamics and directed networks. 
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Controllability of nonlinear systems (including networked systems) has long been an 
important open problem in control theory j5|-0|. Surprisingly strong results were recently 
announced by Wang et al. in Ref. These include a general method to control networks 
of coupled nonlinear dynamical systems, and the full characterization of the minimum 
set of driver nodes. It is proposed therein that any network with one- dimensional node 
dynamics can be fully controlled by the roots of any minimum set of directed trees spanning 
the network, which implies a single driver node in the case of undirected networks (the 
extension to high dimension is immediate Here, we first show that the conclusions in 
Ref. [l[ are generally not valid. Next, we identify and discuss the main flaw that renders the 
proposed control method invalid. We then show that the assertion in Ref. [l[ that networks 
with nonlinear dynamics are more controllable than networks with linear dynamics is both 
not novel and based on wrong reasoning. We also elaborate on the fact that the full control 
of any dynamical system, whether a network or not, concerns much more than merely 
driving it to a target state — it also concerns the stabilization of the target state, which in 
general requires more driver nodes than to reach that state, unless the target state is stable 
[2|. Finally, we comment on the essential literature omitted in Ref. 



(i) Minimum number of driver nodes. 

Arguing that the control of a nonlinear system can be achieved by "linearizing the nonlinear 
system about [a] finite number of local desired states," Ref. [l[ claims that the controllability 
of a network with nonlinear node dynamics can be determined by the structure of so-called 
"local effective networks" (LENs). For one-dimensional nonlinear node dynamics, this leads 
to the conclusion that "any bidirectional network system can be fully controlled by a single 
driver node, regardless of the network topology" [lj]. While these claims are stated for rather 
general systems, networks of n diffusively coupled dynamical units of the following form are 
used as a model system throughout Ref. [lj]: 

n m 

±i = Fi(xi) + ^ a iji H ( x j) ~ H{xi)] + ^ hkU k {t) 

3=1 k=l 
n m 

= Fi{ Xi ) - ^ H ( X j) + J2 b ikMt), (!) 
3=1 k=l 

where x = [xi] nx i and xi stands for the state of node i, which we assume to be one- 
dimensional for convenience of presentation. Here, F^ represents the intrinsic (nonlinear) 
dynamics of node i, H represents the coupling function, A = [ay] rexn is the adjacency 
matrix of the network, L = [£ij] nX n is the Laplacian matrix of the network, and B = [bij] n xm 
is the control matrix (i.e., Uj(t) is the j'th time-dependent control signal). Following Ref. [lj], 
a node i in the network is a driver node if the node is directly controlled, that is, there is 
some k such that bik is nonzero. (Because an irreducible set of s independent input control 
signals can always be assigned one-to-one to s different nodes even if other nodes in the 
network receive a combination of the same inputs, it can be argued that s is more important 
than the number of driver nodes — number of nodes for which "there are input signals on 
them" [l|. But we show below that the results in Ref. [l| do not hold when expressed in 
terms of s either.) 

System ([T]) is controllable if for every given initial state x^ and target state x^ in the 
state space there exists a piecewise continuous input signal u(t) = [uj(t)] mx i such that the 
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system can be steered from x^ to x^ in a finite time. Therefore, the notion of controllability 
is a global property of the system for a given input matrix B. This property has a local 
counterpart, referred to as local controllability. Given a matrix B, the system is locally 
controllable at a given state x^ if there is a neighborhood around this state such that 
every point in this neighborhood can be reached from x^ by a control trajectory in a finite 
time. The authors' identification of the minimum sets of driver nodes required to achieve 
full control of system ([I]) is based on the linearization of the dynamics. However, we show 
that even within the specific class of systems explicitly considered in Ref. the linearized 
dynamics cannot determine the minimum number of driver nodes, let alone their identities. 

Counterexamples to the claims in Ref. [l[ are abundant, even for networks with only two 
nodes. For instance, consider the system 

±i = 2x\ + a 12 (xl - xj) + b\U\(t), 

(2) 

x 2 = x1 + a 2 i(x\ - x\) + b 2 u 2 (t), 

which corresponds to a two-node network in which H{x) = x 2 (see Fig. Ufa)). System is 
controllable if both b\ and b 2 are nonzero (i.e., all nodes are driven). The case of interest is 
the one in which only one of the fej's is nonzero (one driver node). For example, let &2 = 
(similar result holds for b\ = 0). In this case, X2 cannot be decreased from any initial state in 
the state space if 021 G (0, 1), and thus system (J2J) is neither globally nor locally controllable. 
Since the lack of controllability holds for all 021 G (0, 1), system (j2J) is not controllable with a 
single driver node, even though the network has only two nodes. Note that one would indeed 
have reached the opposite conclusion by analyzing the associated LEN as proposed in Ref. 

the linearized dynamics around an arbitrary state = (xi\x 2 ^) T 




2(2 - a 12 )4° 2a 12 4° 
2a 21 x? ) 2(1 - a 21 )xf 




(3) 



and the LEN is a network with edges corresponding to the nonzero entries of matrix J. 
For the range of parameters < a 2 \ < 1 and xf' 7^ 0, for which we have just shown that 
system (J2J) is not even locally controllable, the corresponding LEN is a two-node network 
containing a directed edge from node 1 to node 2 and a self-loop at node 2 (the other edges 
are irrelevant), and is thus predicted by the analysis of Ref. [l[ to be structurally controllable 
through a single control signal input at node 1. This contradiction invalidates the use of 
LEN for the analysis of controllability in networks of nonlinear dynamics. 

In fact, counterexamples are common even when the coupling functions are linear (e.g., 
H(x) = x). Consider a star network with three nodes, with the hub being the only node 
where the control signal is applied (Fig. DJb)). The network dynamics is expressed as 

±1 = Fi(xi) + ai 2 (x 2 - xi) + a 13 (x 3 - x x ) + u^t), 
X2 = xl + a 2 i(xi - x 2 ), (4) 
x 3 = 2x\ + a 3 i(xi - x 3 ), 



< 



which, according to the theory in Ref. [1], is fully controllable for any node dynamics F\ 
for almost all a 2 \ and a 3i such that 021031 > 2 (021 — 031 ) 2 . However, all networks within 
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this range of parameters are in fact not controllable and this follows from Theorems 2.1 
and 2.2 in Ref. [8[. Therefore, the statement that "the minimum number of driver nodes to 
achieve full control of the system is determined by the structural properties of the LENs" 
[l[ is false in general. As further elaborated below, this is so because the controllability 
of the linearized dynamics does not imply (even locally) the controllability of its nonlinear 
counterpart. 

We also comment on the controllability of networks with nonlinear dynamics in gen- 
eral. As shown by Theorem 1 in Ref. [9| and by Theorem 3.1 in Ref. a nonlinear 



dynamical system with n state variables may not be controllable for any given number of 
independent input control signals smaller than n. For example, system §2§ with a single 
input control signal, U\(t) = U2{t), is not controllable for any ai2, 021, 61, and 62 in the 
range [61021 — (2 — 0.12)^2] ■ [&i(l — 021) — ^2^12] > 0, even though both nodes are driven. 
Therefore, this example also rules out the possibility of recasting the Ref. [H's statement 
in terms of the number s of input signals instead of the number of driver nodes, since it 
shows that bidirectional networks with one- dimensional nonlinear node dynamics are not 
guaranteed to be controllable by one input signal, either. On the other hand, the possibility 
of controlling a particular nonlinear system with a small number of inputs is not ruled out. 
The controllability of complex networks with nonlinear dynamics thus depends on both the 
network structure and the form of the dynamics, in sharp contrast to what is stated in 
Ref. 1. 



(ii) Failure of the control method. 

The control method proposed in Ref. [l| is based on linearizing the nonlinear system about a 
finite number of local desired states. The reason this approach fails might appear intriguing, 
since the linearized dynamics approximates the nonlinear system in sufficiently small neigh- 
borhoods of each state. However, when the linearized system is driven between two states 
in one such neighborhood, the resulting trajectory is not necessarily limited to that neigh- 
borhood. This is shown below for linear systems in general, which include as special cases 
the linearized dynamics considered in Ref. [l[. In fact, for such systems, even if the control 
trajectory to go from state x^ to state is strictly local, in general the control trajectory 
to go from state x^ to state x^ will be necessarily nonlocal whenever the number of driver 
nodes is smaller than n. The implication of this is clear: the controllability of the linearized 
dynamics relies on "global" control trajectories and hence this form of linearization cannot 
be reliably used to infer the controllability of the original nonlinear system, even locally. 
To formalize this observation, consider the linear control system in the form 

x = Ax + Bu(t), (5) 

where x = [x±, . . . , x n ] T G M. n is the vector of state variables, u = [u\, . . . , u m ] T is the vector 
of time-dependent control signals, and A = [Aij] nX n and B = [Bij] nxm are constant matrices 
(B plays the role of the input matrix B in system ([1])). 

Definition 1. (Strictly local controllability) - The linear system ([5]) is strictly locally control- 
lable at x^ if for any e > there exists 5 > such that, for any target state x^ satisfying 
ll^ 1 ) _ < 5, the following property holds: there exists a piecewise continuous input 
signal u(t) defined on a finite time interval [to,ti] under which the solution of © satisfies 
x(t ) = x(°\ x(ti) = x^\ and \ \x(t) — x^\\ < e for all t < t < t\. On the other hand, the 
system is not strictly locally controllable at x^ if there exists an e > 0, such that for any 
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5 > there is a target state x^ with \\x^' — x^°'\\ < 5 for which the following property 
holds: any control input that steers the system from x^ ' to x^ 1 ' along a trajectory x(t) for 
to < t < ti with x(t ) = x^ and x(ti) = x^ 1 ' contains at least one point x(r), for t < r < ti, 
such that \\x(t) — x^\\ > e. (Our definition of strictly local controllability is analogous to 
the concept of small-time local controllability often used in the context of nonlinear systems 

Ho 

Figure |2] illustrates this definition. Note that a state of the system can be locally control- 
lable (i.e., it can be driven to any other state in a sufficiently small neighborhood) without 
being strictly locally controllable. A notable exception is when the state is a fixed point of 
the dynamics, as it can be proved that the fixed point of a linear system (which is unique 
and at the origin for nonsingular matrix ^4.) is always strictly locally controllable. On the 
other hand, unless all the nodes are taken as driver nodes, almost every other point in the 
state space is not strictly locally controllable. 

Proposition 1. (Breakdown of strict local controllability) - Consider the linear control sys- 
tem If there is at least one component Xk that is not directly controlled (i.e., Bkj = 
for all j), then the system is not strictly locally controllable almost everywhere. 

Proof. If Akj = for all j, the system is clearly not controllable. If Akj ^ for some j, 
we define the {n — l)-dimensional hyperplane % = {x : (Ax)k = 0}, where we have used 
the notation (Ax)k = Y^j=i AkjXj- The hyperplane % separates the state space W 1 into two 
disjoint regions: "H + = {x : (Ax)k > 0} and %~ = {x : (Ar) fc < 0}. Let the initial state be 
x (o) e -}{+ Choose s > such that 5 e (x (0) ) (lU = 0, where 5 e (x (0) ) denotes the (open) ball 
centered at x^ with radius e. This also implies B £ (x^) n"H~ = 0. For any 6 > 0, choose a 
target state x^ defined by x\ = x[ if i ^ k and x[ = xf^ — 5/2 if i — k. It follows that 
ll^ 1 ) _ a;( )|| = §/2 < 5. In order for the trajectory x(t) with x(to) = x^ to reach x^ in 
finite time (i.e., EUi < oo, such that x(ti) = x^), it is necessary that x(r) G H~ for some 
to < r < ti, since otherwise Xk = (Ax)k would be nonnegative and Xk would not decrease. 
Since B £ (x^) D H~ = by definition of e, we have x(r) ^ £> e (x(°' ) ). See Fig. [3]^a) as an 
illustration. Similarly, we can show that any x G %~ is not strictly locally controllable. 
Since M n = H + U T-L~ U "H and "H is a set of measure zero, the proposition follows. □ 

Figure El^b) illustrates this general result for a network consisting of two nodes with one- 
dimensional intrinsic dynamics. With the exception of the line x± — 0, all the states of this 
system are not strictly locally controllable if node 1 is taken as the only driver node. 

Now, let us return to the specifics of what is done in Ref . [l| • They write the variational 
equation representing the linearization of the nonlinear network dynamics ([T|) around a local 
desired state (Eq. (1) in [l[) in the form 

i = [DF(x (l) ) -L® DH(x®)] ■ £, (6) 

where £ denotes the deviation vector such that £(t = 0) = a^ -* — for the system at the 
local initial state x^ at t — 0. They then go on to write the control problem (Eq. (2) in 
111) in the form 

Z = GZ + Bu(t), (7) 

where G represents the matrix inside the brackets in Eq. OH]). Their approach is illustrated in 
Fig. 1 in Ref. [l[ — a figure strikingly similar to Fig. 1(b) in Ref. [2| — and is then summarized 
as follows: "Once the variational states £ are controlled to approach zero values, the original 
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nonlinear coupled system moves to the local desired state. The next local desired state can 
then be chosen. The final desired state can be reached by repeating this process. Due to the 
strongly structural controllability at each step, the whole controlling process is also strongly 
structurally controllable." 

Before proceeding we note that Eq. ([6]) is incorrect. If the linearization was really around 
the fixed state (as implied in their text), as opposed to the trajectory of x®, then the 
resulting linearized equation would have the extra additive term F(x^) — L® H{x^). This 
term vanishes, in agreement with Eq. flS]) (Eq. (1) in [![), only in the exceptional case when 
x® is a fixed-point solution of the uncontrolled dynamics. In typical cases, when x® is not a 
fixed-point solution, the control of this linear system to reach x® is equivalent to the control 
of a homogeneous linear system away from the origin. For example, when G is invertible, 
the coordinate transformation r\ = £ + [DF(x^) — L ® DH(x^)]~ 1 [F(x^) — L <g> H(x®)] 
brings the problem to the form of Eq. (jSj) but with the target state at [DFfx' 1 ') — L ® 
DH(x^ l ')]~ 1 [F(x^) — L® H(x®)]. The proposed method would then fail because, as shown 
in Proposition 1, the system is not strictly locally controllable at those points and hence the 
linear approximation is violated by the control trajectory (except for the very special cases 
in which all local desired states are fixed points). If instead we accept the authors' form 
of the linearized dynamics (i.e., with no additive term), then £ in Eq. must represent 
the time evolution of the deviation between the trajectories starting at and x^°\ which 
means that the argument of DF and DH should in fact be the time evolution of x® rather 
than itself. In Ref. [l[ it is claimed (as a crucial condition within their approach) that 
matrix G is a constant, which, in the best case scenario, is only an approximation since G 
varies along the trajectory of x®. This can be criticized in its own right, but let us assume 
for the moment that this approximation is acceptable. 

One might then be tempted to reason that, because the approach in Ref. [lj is based on 
a linearization around the origin (with respect to £), the breakdown of strict local control- 
lability would not come into play. However, unless x® is a fixed point of the uncontrolled 
dynamics, £ being controlled to zero in a finite time T does not mean that the control tra- 
jectory will reach the desired state x® — this only implies that the control trajectory will 
meet the uncontrolled trajectory of x"> at time T. While one can reduce the convergence 
time T arbitrarily and hence bring the meeting point to be arbitrarily close to the desired 
state x® (at least within the framework of unconstrained control), in general this necessar- 
ily causes the control trajectory to go arbitrarily far away from the origin (except for the 
unlikely situation in which is zero for the components of all nodes that are not directly 
controlled, i.e., when they are already equal to those of the desired state). Therefore, the 
trajectory is necessarily nonlocal in all nontrivial cases and the linearization becomes an 
invalid approximation of the nonlinear dynamics. This can be formalized as follows. 

Proposition 2. (Time lengthening versus trajectory diversion to reach the origin) - Consider 
the linear control system fl5]). Suppose that there is at least one component Xk that is not 
directly controlled (i.e., Bkj = for all j). Let x^ = (x[ , . . . , x„ ) be the initial state and 
x^ = (0, . . . , 0) be the target state. Let x(t) be a control trajectory such that x(0) = x^ ' 
and x(T) = x^\ and let e = max <t<T Then T > \xl \-Akj\s)- 

Proof. Along any control trajectory x(t) with ||x(t)|| < e, we have |i/c(t)| = | ^2jAkjXj\ < 
Y^a |^fcj|rnaxj \xj\ < |*4;y|£. Since Xk(0) = xf^ and Xk(T) = 0, it follows that \x^\ < 
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In the case of a network with one- dimensional node dynamics, having at least one 
component Xk that is not directly controlled corresponds to having a number of driver 
nodes smaller than n (the case of a single driver node is, of course, the most extreme of 
all such cases). Therefore, as anticipated above, it follows immediately from Proposition 2 
that in general any attempt to make T small will lead s = maxo<j<T | \x(t) \ | to be very 
large. When e is large, the problem becomes nonlocal and linear approximation fails, as in 
the case of non-strictly local controllability established in Proposition 1. 

(Hi) Failure of variants of the control method. 

The presentation in Ref. [l[ involves a number of contradictions. These include the inconsis- 
tent definition and use of £ and the use of x l (represented by x^ here) to mean incompatibly 
different objects in the text, equations and figures (sometimes different things in different 
parts of the same equation). It is thus instructive to clarify that alternative interpretations 
of the authors' purported method are also destined to failure. 

Specifically, since x^ cannot in general be brought to a given point x"' with the lin- 
earization procedure proposed in Ref. [lj, one may ask whether steering the trajectory of 
x^ towards a different point of the trajectory of would be any more promising. Needless 
to say, if control can only bring the trajectory of x^ to meet the trajectory of x^ at a point 
in the forward part of the trajectory with respect to x"' itself, then the system cannot be 
controlled in a direction that goes against the flow (except for very special cases). The case 
in which the meeting point would be in the backward part of the trajectory of is po- 
tentially more promising because then the system could evolve spontaneously to the desired 
state x^ l \ However, we now show that, unless the desired state is in a small neighborhood 
of a fixed point of the original (nonlinear) dynamics, this is generally not possible for most 
initial states without violating the linear approximation. 

Proposition 3. (Local drift of uncontrolled variables) - Consider the control system in the 
form x = F{x) + Bu(t), where J 7 can be nonlinear, and assume that Xk is not directly 
controlled (i.e., Bkj = for all j). Let be a state such that J c (x)k > for all i in a 
neighborhood B e (x^) for some e > 0. A state x^ E B £ (x®) cannot be reached from any 
initial state x^ E B e (x®) with x£ > x k through a control trajectory inside B e (x®). 

Proof. By assumption, ik = J r (x)k > for all x E B e (x®), regardless of the input u(t). 
It follows that any trajectory x(t) inside B £ (x^) with a; (to) = necessarily satisfies 
x k {t) > xf > x[ l) for all t > to. □ 

This proposition implies that, in general, for most initial states x^ in the neighborhood 
of x^ l \ one cannot reach any point a^ 1 -* in the backward trajectory of x® by means of a local 
control trajectory. This property underlies the breakdown of strictly local controllability in 
both linear and nonlinear systems. 

(iv) Linear versus nonlinear systems. 

It is claimed in Ref. [H that "network systems with nonlinear dynamics are more control- 
lable than with linear dynamics." This conclusion is based on the authors' interpretation 
that for networks with linear dynamics "there are limitation [s] on the network controllabil- 
ity, depending on the particular network structure," while networks with (one-dimensional) 
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nonlinear dynamics can be "fully controllable by one controller at any node, regardless of 
the network structure." We have already demonstrated that their conclusion about net- 
works with nonlinear dynamics is incorrect. We now show that their interpretation of the 
controllability of networks with linear dynamics is also incorrect. 

The authors' criterion is entirely based on the minimum number of driver nodes to control 
the network. According to this criterion, however, networks with nonlinear dynamics cannot 
be more controllable than those with linear dynamics for the class of systems considered in 
Ref. This is so simply because connected bidirectional networks with one-dimensional 
linear dynamics are always controllable by a single driver node for almost all weights of 
the edges whenever the intrinsic node dynamics is non-degenerate, as expected in general. 
This has been clearly demonstrated in Ref. where the case of directed networks was 



also properly considered, and can also be anticipated from Ref. [12[. Within the model 
system ([1]), considered in Ref. [H, the linear case corresponds to taking Fj(xj) = CjXj and 
H(xi) = Xi, where q is a constant, and the non-degeneracy condition corresponds to at 
most one q being zero (this assures that pi = Ci — Ylj^i a ij are a ^ distinct for almost all 
choices of edge weights). The results of Ref. [3| are in fact strikingly similar to those of 
Ref. [H, except for the fact that the former are prior and correct while the latter are not. 
Incidentally, despite all this, Ref. [3[ is not referenced in Ref. 

Having shown that, according to the criterion of Ref. networks with nonlinear dy- 
namics could not be "significantly more controllable" than networks with linear dynamics, 
we clarify that for the purpose of this comparison their criterion itself is misleading. Control 
concerns not only bringing the network to a specific state, but also keeping it at that state. 
Keeping the network at a state other than a stable state generally requires direct control of 
all nodes (a known fact also discussed in Ref. |H). This is illustrated in Fig. E^b), where the 
system can be brought from x^ to x^ by controlling only one of the two nodes but the 
direct control of both nodes is required to keep the system at x^ . This is a very general 
and widely appreciated property in control theory. 

While the linearization procedure and the concept of LEN proposed in Ref. [l| do not 
contribute to our understanding of the controllability of nonlinear systems, other valid in- 
roads into this area had already been announced and are ignored in Wang et a/.'s account 
of the literature. An obvious precedent ignored in Ref. [l| is Ref. 0, where some of us de- 
veloped a general method to control networks with nonlinear dynamics under rather general 
conditions (including arbitrary constraints on the control interventions, which are expected 
to be crucial in realistic situations). The method in Ref. [2[ drives the system to the basin 
of attraction of a desired stable state (from which the uncontrolled system can evolve spon- 
taneously to this desired state), and hence does not depend on continuous control to keep 
the system in that state. The existence of multiple stable states, including desirable ones, 
is possible due to the nonlinear nature of the dynamics (nonsingular linear systems have at 
most one stable state) and is known to be common in real complex networks of interest, 
such as food webs, metabolic networks, and power grids. Within this framework, it was 
already noted in Ref. [2[ that "The possibility of directing a complex network to a prede- 
fined dynamical state [...] can be achieved under rather general conditions by systematically 
designing compensatory perturbations that [...] take advantage of the full basin of attraction 
of the desired state, thus capitalizing on (rather than being obstructed by) the nonlinear 
nature of the dynamics." This is one criterion (accounting for the fact that the system has 
to be kept at the desired state) in which networks with nonlinear dynamics are indeed more 
"controllable" than networks with linear dynamics p|. Another angle from which the same 
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conclusion has been reached is when the control itself is allowed to be nonlinear, owing to 
independent directions of controlled motion (in addition to those from the linearized control 



dynamics) generated by the Lie algebra of the control vector fields [13H16 



The algorithm presented in Ref. 0] to design control interventions is based on iterative 
applications of the variational equation to identify state perturbations that bring the closest 
approach point of the trajectory incrementally closer to the desired state. Equation (1) 
in Ref. on which their method is built, is an ill-conceived implementation of the same 
variational equation to a situation where, as explained above, it cannot be used. Therefore, 
the whole framework in Ref. [l[ has major problems that go well beyond the missing 
literature and which render their results unusable. A main reason iterative linearization 
can be used in Ref. J2['s formulation (but not in Ref. (H's) is that to reach a target basin of 
attraction all it is needed is that any point of the basin be accessible from the initial state. 
The identification of such points is a challenging problem due to the constraints and the 
high dimensionality typical of large real networks, and this is a central question addressed 
in Ref. Q. 



(v) Final remarks. 

In summary, the results in the paper of W.-X. Wang, Y.-C. Lai, J. Ren, B. Li & C. Grebogi [l[ 
are invalid. While we have focused on networks of one-dimensional dynamical systems, 
which form the bulk part of their work, it is immediate that their generalizations to high- 
dimensional node dynamics and networks with directed edges are also invalid. Unfortunately, 
the control method proposed in Ref. [l| cannot be easily fixed because, as shown here, 
the results of their linearization-based method and local effective networks fail both as a 
necessary and as a sufficient condition for the controllability of the actual nonlinear dynamics 

(Fig. m). 

Note that most of the analysis presented in this comment (in the form of simple but 
rigorouspropositions) is expressed in a much stronger form than needed to point out the flaws 
in Ref. For example, our analysis reveals comprehensive limits for the controllability of 
networks with nonlinear dynamics based on the proposed local linearization for any number 
of driver nodes smaller than n (the case of a single driver node is, of course, the worst of 
all such cases). In particular, the present arguments show that, based on local linearization, 
even for n — 1 driver nodes such networks would be far less controllable than claimed in 
Ref. [l[ for a single driver node. 

Leaving aside the technical problems, the concept of controllability that is the focus of 
Ref. PJ concerns the ability of steering the system from any given state to all other states 
in the state space. In practice, however, this notion is often of limited usefulness since what 
is often needed (in both engineered and natural networks) is control within specific regions 
of the state space. In addition, far more important than the ability to move the system 
between any two states in the state space is the ability to account for constraints on the 
controls that are unavoidable in realistic situations. These constraints often take the form 
of limiting the direction and magnitude of the control at certain nodes, while forbidding 
the direct control of others. Positive advances along these lines can and have been made 
[1, [l8[ . We hope the limitations of the approach commented upon here will not discourage 
others from contributing to this important area of research. 
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FIG. 1: Counterexamples to the assertion that the controllability of a network with nonlinear 
dynamics can be predicted from the linearization procedure proposed in Ref. The solid lines 
represent edges in the network, while dashed self-loops represent additional edges in the LENs 
derived from Ref. (a, b) Examples of networks whose linearized dynamics are controllable while 
their nonlinear counterparts are not: (a) bidirectional network with two nodes corresponding to the 
dynamics in Eq. ([2]) and (b) bidirectional network with three nodes corresponding to the dynamics 
in Eq. @. (c) Example of a network whose nonlinear dynamics is controllable while its linearized 
counterpart is not: star three- node network endowed with the dynamics x\ = u\(t), x<i = 021^1 , 
%3 = a 3l x ii where 021 7^0 and 031 7^ 0. The conclusion that this system is controllable follows from 
Theorem 2.1 in Ref. [17]]. Therefore, the "controllability" determined by the linearization-based 
method of Ref. [l| fails both as a necessary and as a sufficient condition for the controllability of 
the actual nonlinear dynamics. 



(a) ...... (b) 




FIG. 2: Schematic representation of strictly local controllability, (a) Illustration of a strictly locally 
controllable state x^: nearby states can be reached from this state through a control trajectory 
that is limited to a small neighborhood, (b) Illustration of a locally controllable state that is 
not strictly locally controllable: there are nearby states that can only be reached through control 
trajectories that are not limited to a small neighborhood. 
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FIG. 3: Illustration that almost all states of linear systems are not strictly locally controllable, (a) 
Schematic representation of this property as formalized in the proof of Proposition 1. (b) Example 
two-dimensional system x\ = x% + ui(t), X2 = xy. any point not on the line x\ = is not strictly 
locally controllable. The background arrows indicate the vector field in the absence of control. 
The curve represents a typical control trajectory to bring state x^ to state x^\ illustrating that 
whenever these states are in the x\ > half of the plane and £3 < £3 , the control trajectory has 
to visit the x\ < half of the plane. Note that this is so because not all the variables are directly 
controlled, which is tantamount to not controlling all the nodes in a network. 
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